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Abstract 

The percolation of strings gives a good description of the RHIC 
experimental data on the elliptic flow, V2 and predicted a rise on the 
integrated V2 of the order of 25% at LHC such as it has been experi- 
mentally obtained. We show that the dependence of V2 on pj- for RHIC 
and LHC energies is approximately the same as it has been observed, 
for all the centralities. We show the results for different particles and 
the dependence of V2 on the centralities and rapidity. Our results are 
compatible with an small value of the ratio ij/s in the whole energy 
range such as it was expected in the percolation framework. 

1 Introduction 

A major breakthrough was the discovery at RHIC experiments of a large 
elliptic flow V2 [1-9]. The observed anisotropic flow can be understood only 
if the particles measured in the flnal state depend not only on the physical 
conditions realized locally at their production point but also on the global 
geometry of the event. In a relativistic local theory, this non-local informa- 
tion can only emerge as a collective effect, requiring interactions between 
the relevant degrees of freedom, localized at different points of the collision 
region. These anisotropic flow is particularly unambiguous and convincing 
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manifestation of collective dynamics [10]. The large elliptic flow V2 can be 
qualitatively explained as follows. In a collision at high energy the specta- 
tor nucleons are fast enough to move away leaving behind at mid-rapidity 
and almond shaped azimuthally asymmetric overlap region filled with the 
QCD matter. This spatial asymmetry implies unequal pressure gradients 
in the transverse plane, with a larger density gradient perpendicular to the 
reaction plane (in- plane). As a consequence of the subsequent multiple inter- 
actions between the degrees of freedom involved, the spatial asymmetry leads 
to anisotropy in the momentum space [4,11,12]. The final particle transverse 
momentum is more likely to be in-plane that in the out-plane, with V2 > 
as predicted [5] . 

The general idea has been realized in various mechanisms for the source of 
elliptic flow. A convenient and successful way to describe the flow anisotropy 
is achieved in the hydrodynamical approach [6,8], taking into account the 
asymmetric shape of the nuclei overlap in the transverse plane at values 
of the impact parameter b different from zero. The description of course 
assumes collective effects to be responsible of the flow. The microscopic 
framework of the hydrodynamics approach, as well as the understanding of 
the required early thermalization remain open questions. The percolation of 
strings [13,14] can provide the microscopic picture and the answer of these 
questions. 

In a heavy ion collision at high energy, strings stretched between the par- 
tons of the projectile and the target are formed. Each string in the transverse 
space looks like disk of area Trrg, with ro — 0.2 — 0.3 fm due to the confine- 
ment. As the energy and/or the size of colliding objects grows the number 
of strings increases and the strings overlap forming clusters, whose behavior 
is determined by the color field inside which is given by the vectorial sum of 
the color field of the individual string. In this way, the formation of clusters 
can be seen as interactions of the partons of the individual strings with the 
corresponding color rearrangement. On the other hand, a local temperature 
can be defined associated to the fragmentation of the string via the Schwinger 
mechanism [15,16], and in the same way a temperature can be associated to 
the cluster formed from many individual strings. As far as this cluster covers 
most of the total interaction area, (what happens if the string density is over 
the critical percolation threshold) this temperature becomes a global tem- 
perature determined by the string density. This fact, together with the good 
description of the elliptic fiow of RHIC data makes the percolation of strings 
a good candidate to be the initial state for hydrodynamical description. 
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The elliptic flow was studied in the framework of percolation in references [17- 
20] leading to the developments of analytical formulaes for f 2 including the 
dependence on transverse momentum, which were successfully compared with 
the RHIC experimental data [17-18]. In reference [20] were given arguments 
in favor of the methods used previously. 

In this paper we present our results for LHC energies. Previously, we 
already anticipated a rise of around 25% [19] for the integrated elliptic flow as 
it has been obtained experimentally [21]. In the evaluation , the main source 
of uncertainties is the string density due to the uncertainties in the number 
of strings and its dependence on the energy and centrahty. To avoid these 
uncertainties we determine the string density in Pb — Pb at LHC at different 
centralities comparing the corresponding formula in percolation with the data 
of the ALICE collaboration [22] . Given the string density, the elliptic flow is 
obtained using our analytical formula. 

We find that V2 for pr < 0.8 GeV/c is very similar at LHC and RHIC 
energies, being slightly larger for px > 0.8 GeV/c. Concerning the rapidity 
distribution there is no longer a triangle shape and a flatten distribution is 
obtained for \y\ < 4. Also there is not limiting fragmentation in the same 
way that in percolation there is not limiting fragmentation in We obtain 
a very reasonable overall agreement with the experimental data, compatible 
with the existence of a partonic fluid with a low ratio between the shear 
viscosity and the entropy density in the whole energy range from RHIC to 
LHC. This low ratio for this energy range was expected in the percolation 
framework [23-24]. 

The plan of the paper is as follows. In the next section we present a brief 
introduction to percolation of strings. In section 3 we describe the obtention 
of the analytical formulas for the elliptic flow. In section 4 we discuss our 
results comparing with LHC and RHIC experimental data. Finally, in section 
5 we present our conclusions. 

2 The string percolation model 

In the collision of two nuclei at high energy, the multiparticle production 
is described in terms of color string stretched between the partons of the 
projectile and the target. These strings decay into new ones by qq or qq — qq 
pair production and subsequently hadronizc to produce hadrons. Due to 
the confinement, the color of these strings is confined to a small area in 
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transverse space 5*1 = nr^ with tq ~ 0.2 — 0.3fm. This value corresponds 
to the correlation length of the QCD vacuum [25]. With increasing energy 
and/or atomic number of the colliding particles, the number of strings grows 
and they start to overlap forming clusters. At a certain critical density, a 
macroscopical cluster appears, which marks the percolation phase transition. 
This value corresponds to the value Pc = l-18~l-5 (depending on the type of 
the profile functions of the nucleus employed, homogeneous or Wood-Saxon) 
where the variable p = NsSi/Sa, and Sa 



Sa = 2Rl[cos-\P) - Py/T^] (1) 

with 

is the overlapping area of the nucleus. In order to describe the behavior of the 
clusters formed by several overlapping strings we must introduce some dy- 
namics. We assume that a cluster of n strings behaves as a single string with 
an energy-momentum that corresponds to the sum of the energy-momentum 
of the individual strings with a higher color field, corresponding to the vec- 
torial sum of the color field of each individual string. In this way [14] [26] we 
can compute the mean multiplicity /in and the mean transverse momentum 
squared < of the particles produced by a cluster: 



^ V ~s^'^^' ^ "^"^ ^^^^^ 

where pi and < >i are the mean multiplicity and mean pf. of particles 
produced by a single string and Sn is the total area occupied by the n disks, 
which can be different for a different configurations even if the clusters have 
the same number of strings. Note that if the strings just touch each other, 
Sn = nSi, and the strings act independently of each other. In contrast, if 
they fully overlap Sn — Si, and it is obtained the largest suppression of the 
multiplicity and the largest increase of the transverse momentum. In the 
limit of high density, one obtains [14] [26] 

Si P I ,.s 

and equations (3) transforms into analytical ones 

2 

p = N,F{p)p„ < pI >i= (5) 
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In the mid-rapidity region, the number of strings is proportional to the num- 
ber of colUsions Ncoii ~ and grows hke s^^. Therefore p ~ N'^^ , and 
IX ~ Na and ^ s^. In other words, the muhiphcity per participant is al- 
most independent of Na (the only dependence on Na arises from the factor 
Vl — e~P in eq. (5)). Outside the midrapidity, Ng is proportional to Na 
instead of N^^. Therefore, there is an additional suppression factor N]^^ 
compared to central rapidity. 

Concerning the transverse momentum distribution, one needs the distri- 
bution g{x,pT) for each cluster and the mean squared transverse momen- 
tum of the clusters, W{x). For g{x,pT) we assume the Schwinger formula 
g{x,p^) — exp{—p^x). For the weight function, we assume the gamma dis- 
tribution 

^<-) = W^"' 

which is the simplest distribution among the distributions stable under the 
size transformations [27,28]. In (6) 7 and k are 

l_<x^>-<x>'^ _ k , , 

k <a;>2 <x> 

I is proportional to the width of the distribution, depending on the string 
density p. At small density, there is not overlapping of strings and all the 
strings produce particles with the same mean transverse momentum. The 
width is zero and A; — > 00. When p increases there is some overlapping of 
strings, different clusters are formed and k decreases. The minimum of k 
will be reached when there are more different clusters. Above this point, the 
different clusters start to join and k increases. In the limit p ^ 00, there is 
only one cluster formed by all the produced strings and A; — > 00. Therefore 
the transverse momentum distribution /(pt, u) of the particle i is given by 
[28] 



dN 

fiPt,y) = , 2j = d{x)W{x)g{pT,x) 



(8) 



where the dependence of A; or p is determined by comparing (8) with the 
experimental data. The formula (8) is valid for all the energy, centralities 
and pt < 4: — 5 GeV/c. In the case of baryons we must introduce minor 
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changes [29]. The equation (8) is the main ingredient for the evaluation of 
the eUiptic flow. 



3 Elliptic flow 



The cluster formed by the strings has generally an asymmetric form in the 
transverse plane and acquires dimensions comparable to the nuclear overlap. 
This azimuthal asymmetry is at the origin of the elliptic flow in percolation. 
In fact, the partons emitted at some point inside the cluster have to pass 
a certain length through the strong color field of the cluster before they 
appear outside. It is natural to assume that the energy loss by the parton is 
proportional to the length and therefore the transverse momentum pt of an 
observed particle will depend on the path length travelled and so different 
for different direction of emission [19]. 

In this way the inclusive cross section ^^i^y^^ = fiPri PipiV) ob- 
tained from eq (8). doing the change p ^ p^, where 




(9) 



where R and R^p are shown in fig 1. Notice that 




(10) 



Expanding /(pr, P^iV) around p or i? we have 



f{PT, p<f, y) ^ -/(pt, p, y)[i H ^ 




(11) 



Notice that due to (10) we have 




(12) 



and from (11), finally we obtain 




9 ^ I 

V2{PT,y) = - dipcos{2ip)[l 



R ' ' 2F{pf > (1 + F{p)pl/ < pI >i) 



(13) 
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Figure 1: Azimuthal dependence of R. 



having present that 



sin((y9 — a) 



(14) 



Ra simp 
a = sin~^(/3sin(/9) 



(15) 



for fixed ^/s and p, in the hmit cases of 6 — ?> and b — 2Ra we have: 

(i) 6 — )■ (i.e /3 — )■ or Na A) which imphes a — )■ 0, and V2{pt, u) 0. 
{ii) b 2Ra (i.e (3 ^ 1 or Na — )■ 0) which imphes a = (p and W2(Pt, I/) 

if we look now to dependence of V2, we see that ^2 — > as — > and 

V2 — )■ constant as p| — > oo. 

We observe that at low pt, say pt < 1, the dependence on p of (13) is 

given by {e~P — F{pY) /2F{p) which remain approximately constant between 

p = 2 and 5 which corresponds to central Pb — Pb collisions at RHIC and 

LHC energies, therefore we expect the same V2 at low pr at RHIC and LHC 

energies. 

We perform next the integration in weighted by ir^l^ ^o obtain 



Again, we have ^ as Na A and also ^2 ^ as Na 0. As the 




(16) 
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bracket in (16) is in modulus a growing function of p, V2 at fixed rapidity is 
a growing function of energy and 

In order to obtain formula (13) we have retained only the first term of the 
expansion on R^. This is a reasonable approximation at low pt- For higher 
Pt we expect corrections of the order of 10 — 15%. 



4 Results and discussion 

The first formula in (5) can be written as 

rr- 7 r— TTTq 1 du 

^(l_exp(-p))p = -°— — (17) 

/^i aV 

As we know ro = .2 fm and /^i = .8 obtained previously by comparing eq. 
(17) with pp data, we can compute the values of p from eq. (17). This 
allows us to eliminate the main source of uncertainties which comes from the 
value of the number of strings Ng and its dependence on the energy and on 
the centrality. We know that A^^ ~ s^^ and A^^ ~ ^-iid we obtained 

2 A ~ 2/7 from energy-momentum conservation arguments [34], but in this 
paper, we determine the string density directly from the experimental data 
on ^ to eliminate any uncertainty. Prom the centrality dependence of the 
experimental data [22] at a/s — 2.76 TeV for Pb — Pb coUisions and at 
^/s = 200 GeV for Au — Au collisions we obtained the dependence of p on 
the centralities at these two energies as it is shown in Fig. (2). In Fig. (3) we 
show the obtained centrality dependence oi with the p- values together 

with the experimental data [22] . 

Prom the eq (13), we compute the transverse momentum dependence of 
V2 at different centralities. The results of f 2 at different centralities are shown 
in figures (4) to (7), for Ph - Ph a.i = 2.76 TeV compared to the RHIC 
results for Au - Au with ALICE and STAR data [30]. We obtain at all 
centralities similar values of V2 for < -8 GeV/c and for pt larger the LHC 
values are shghtly larger (10%) than RHIC values. The difference is of the 
same size than the experimental ones. Very similar results were obtained in 
the percolation framework [20] using different approximations. In this case, 
the values found are the same at RHIC and LHC in the whole pr range. 

At = 5.5 TeV, from the formula (13) we expect that V2 for low px 
would be slightly smaller than at ^/s — 2.76 and almost equal at higher pt- 
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Figure 2: String density dependence with the number of participants, dashed 
and sohd hues correspond to i/s = 2.76 TeV and i/s = 200 GeV energies 
respectively. 




Figure 3: hues in red and blue correspond to the dependence for the a/s = 
2.76 TeV and ^/s = 200 GeV energies, errorbars in green and pink correspond 
to ALICE data [22] and PHENIX data [31] respectively. 
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Figure 4: Stars in red and blue squares correspond to our predictions for 
^/s = 200 GeV and ^/s = 2.76 TeV energies, and errorbars in green and pink 
are the respective data from references [22] [30] for centrality 10 — 20%. 
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Figure 5: Stars in red and blue squares correspond to our predictions for 
^/s = 200 GeV and ^/s = 2.76 TeV energies, and errorbars in green and pink 
are the respective data from references [22] [30] for centrality 20 — 30%. 
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Figure 6: Stars in red and blue squares correspond to our predictions for 
^/s = 200 GeV and ^/s = 2.76 TeV energies, and errorbars in green and pink 
are the respective data from references [22] [30] for centrality 30 — 40%. 
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Figure 7: Stars in red and blue squares correspond to our predictions for 
^/s = 200 GeV and ^/s = 2.76 TeV energies, and errorbars in green and pink 
are the respective data from references [22] [30] for centrality 40 — 50%. 
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Figure 8: Results for elliptic flow dependence with pseudorapidity in Pb — Pb 
collisions at ^/s = 2.76 TeV with Npart = 211. 



In Fig. (8) and (9) we plot the rapidity dependence of V2 for central 
Pb — Pb collisions at = 2.76 TeV and for Au — Au collisions at y/s = 200 
GeV at Npat = 211 respectively, obtained from eq(16). In Fig. (9), we 
include the corresponding experimental data [32]. The integrated V2 obtained 
is around 25% higher at LHC than at RHIC energies in agreement with the 
experimental data [21]. We observe a breaking of limiting fragmentation as 
it is expected for percolation in ^ [33]. 

In Fig. (10), we compare the values of V2 at different centralities with 
the experimental data, obtaining a good agreement. In fig (11) we show the 
results at = 2.76 TeV for vr, k, and p at a centrality of 40—50% for Pb—Pb 
collisions. In figure (12) we check the constituent quark scaling plotting the 
ratio between the elliptic flow and the number of constituents as a function 
of the ratio between the kinetic energy and the number of constituents at 
■\/s = 2.76 TeV. In this case, we observe an approximate scaling at low values 
of the kinetic energy. For high values, we have deviations of the scaling at 
y/s = 2.76 TeV of the order of 20%. In our evaluations, we have neglected 
terms that for p-p > 1 GeV/c can give contributions of the same order. 

The similar values of f 2 found for LHC and RHIC energies points out the 
existence of a liquid of very low value of the ratio between the shear viscosity 
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Figure 9: Triangles show the corresponding resuhs for elhptic flow depen- 
dence with pseudorapidity in Au — Au colhsions at ^/s = 200 GeV with 
Npart = 211, compare to data from PHOBOS [32]. 
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Figure 11: Red doted line, green solid line and blue dashed line are correspond 
to the proton, kaon, and pion predictions for central Pb — Pb collisions at 
v/i = 2.76 TeV. 



1] and the entropy, s. This low ratio is expected in percolation approach [23]. 
In fact, the ratio is given by 

^ - ^ (18) 

s hnatr 

where T is the temperature, n the number of effective sources per unit volume 
and atr is the transport cross section. According to the second part of the 
formula (5), the transverse correlation length is ^jF{p) and therefore 

atr = 7^rlF{p) = S^F{p) (19) 

On the other hand, the number of sources is given by the ratio between 
the area covered by strings and the area corresponding to the transverse 
correlation length 



and therefore 



with L ~ 1 fm. 



14 



0.14 



0.12 



0.08- 




0.06- 



0.04 



0.02 



KET/n„ 



Figure 12: Elliptic flow scaling with the number of quark contituents for 
central Pb — Pb colhsions at ^/s = 2.76 TeV. Symbols (+) in red ,(*) in blue 
and (x) in green are used for protons kaons and pions respectively. 
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Finally, we can relate the temperature to the tension of the cluster [16] 
resulting 

T = {<pI >i /2F(p))V2 (22) 
Introducing (19), (20), (21) and (22) into (18) we obtain 
77 ^ <pt>iL p^'^ 

s 5V2 (l-e-^)5/4 ^ > 

where we have approximated (< >i)^^^ — < Pt >i- For < Pt >i= 200 
MeV/c, close to the threshold of percolation p ^ 1.2 — 1.5, we obtain rj/s c::^ 
0.2 and for p ~ 3 — 5 we obtain rj/s ~ .25. The ratio grows slowly, as 
pi/4 remaining low between RHIC and LHC energies. Let us remark, that 
this low ratio is due to the small transverse correlation length, decreasing as 
the energy increases compensating the increasing of the effective number of 
sources in such a way that the mean path length X — 1/ na — remains 
almost constant. 



5 Conclusions 

The close analytical formulas obtained in the percolation framework from 
the approximated transverse momentum distribution allow us to compute 
the elliptic flow once the string density is known, which is fixed from the 
experimental data at LHC and RHIC energy. The transverse momentum 
dependence of V2 and its values are very similar at RHIC and LHC energies 
for all centralities, in agreement with the experimental data. The integrated 
elliptic fiow at mid rapidity is around 25% higher at a/s — 2.76 TeV than at 
^/s — 200 GeV and in the rapidity dependence of the elliptic flow there is no 
longer a triangle shape shown at RHIC. In the percolation framework there 
is not limiting fragmentation for both V2 and dN/dy. The whole picture 
is consistent with the formation of a fluid with a low ratio between shear 
viscosity and entropy in the range of energies of RHIC and LHC such as 
it is expected in percolation. The percolation framework provide us with a 
microscopic partonic picture which explains the early thermalization and the 
large interaction required by the hydrodynamical approach. 
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